Pairwise Summation Approximation of Casimir energy from first principles 
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We obtain the Pairwise Summation Approximation (PSA) of the Casimir energy from first principles in the 
soft dielectric and soft diamagnetic limit, this analysis let us find that the PSA is an asymptotic approximation of 
the Casimir energy valid for large distances between the objects. We also obtain the PSA for the electromagnetic 
(EM) coupling part of the Casimir energy, so we are able to complete the PSA limit for the first time for the 
complete electromagnetic field. 
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I. INTRODUCTION 

Since 1948, when Casimir introduced the energy that got 
his name IH, calculation formulas have been looked for Many 
analytical and numerical methods have been proposed, such 
as the zeta function technique, the heat kernel method, semi- 
classical methods or Green function (local) methods just to 
mention a few of them |2]. However, exact results have been 
obtained only for some simple geometries. 

The asymptotic analysis of Casimir energies has also a long 
history. In fact, Casimir himself and Polder in year 1948 gave 
the first asymptotic formula for the Casimir energy between 
two electrically neutral bodies in terms of their electric in- 
duced dipoles jst). Some time later, a generalization of that 
formula, known as the Pairwise Summation Approximation 
(PSA) was derived |4] for electric media. The main assum- 
ption is a linear superposition of the Casimir - Polder interac- 
tions between the induced polarizabilities of each element of 
volume body. Then, the PSA energy is expressed as an inte- 
gral over the two object's volumes and it is proportional to the 
objects polarizabilities. The formula has been recently reob- 
tained by R. Golestanian in and by K. A. Milton et. al. in 
[6] in the soft dielectric limit. 

Besides, a new asymptotic method for calculating Casimir 
energies in term of the induced multipoles of the interacting 
bodies has been proposed in 1 7 ] and [ 8] . This formula provides 
a procedure for the calculation of the Casimir energy between 

arbitrary shaped compact objects ^9^. 

Our goal is the asymptotic calculation of the electromagne- 
tic Casimir energy and a presentation of a systematic asym- 
ptotic expansion procedure for the integral dipole formula on 
higher orders. We will reobtain the classical results of Casimir 
and Polder |3], and of Feinberg-Sucher in jlOll in the rarified 
or soft dielectric and soft diamagnetic limit. 

The result we obtain here is a generalization of Milton et. 
al. one [6] but assuming that the bodies are soft diamagne- 
tic as well as soft dielectric. For this purpose, we will use the 
multi-scattering expansion of the Casimir energy formalism 
given in 1 7]. To our knownledge, it is the first time this forma- 
lism is used to derive the PSA. We obtain that the PSA is the 
first order of a perturbation expansion in the difference bet- 
ween the electric and magnetic permeability constants of the 
objects respect the electric and magnetic permeability cons- 
tants of the medium were the objects are placed. The interest 
of this result is that we obtain this pairwise summation formu- 



la for the complete electromagnetic Casimir energy which is 
derived as an asymptotic limit of an exact and free of divergen- 
ces formula. That means that now we can establish the range 
of validity for that approximation. In fact, we will justify why 
this formula is valid in the far distance objects limit. This de- 
rivation also gives us a perturbative procedure for corrections 
of this approximation and posterior expansions to more than 
two objects or finite temperature cases. 

We will follow this plan for the article: Using the soft die- 
lectric and soft diamagnetic approximation, in Sect.|ll]we will 
obtain the PSA of Casimir energy for the zero temperature ca- 
se starting from the exact Casimir energy formula given in 
We will also study the far distance limit to reobtain the asym- 
ptotic Casimir energies given in |3] and [10]. In Sect. Hill we 
will obtain PSA Casimir energy formulas for any temperature, 
studying the high and low temperature limits. In Sect. llVlwe 
will study the PSA Casimir energy for the three bodies system 
and for the general N bodies system. We will obtain a kind of 
superposition principle of Casimir energy in the PSA approxi- 
mation. Finally, in Sect. [V]we will study the first perturbation 
energy term to the PSA approximation and we will discuss 
about the nature of the PSA limit of the Casimir energy. 



II. DILUTED LIMIT AT ZERO TEMPERATURE 

Our objective is the calculation of the complete electromag- 
netic Casimir energy between two bodies in the soft dielectric 
limit. For this purpose, we will use the Casimir energy formu- 
la between two compact bodies at temperature T — 0, given 
in H as 
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(11. 1) 



Where k is the frequency, I is the identity matrix, N is the 
matrix N = ¥111121211121. Here, is the T scattering matrix 
of the a - body under the electromagnetic field, and Vap is 
the propagation matrix of the electromagnetic field from a 
- body to [3 - body. We will use the position representation 
instead the multipole representation used in |7], so we can 
identify Uap ~ Goap, where Goap is the free dyadic Green 
function. Taking into account that log(|A|) = Tr(log(A)) and 



2 



that log(l 



transform Eq. ( III. Il l into 
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Equation ( III.2I I is an asymptotic expansion of Eq. (III. Il l in the 
Pmax < oo case, which means that our calculus will be va- 
lid in the large bodies distance limit. The T operator is related 
with the potential V (we will see what V is later) by the Lipp- 
mann - Schwinger equation, that can be written as li 1 li1 



T = (1 - VGo) ^ V. 



(11.3) 



Applying a Born expansion to Eq. (|II.3l l, we can obtain an 
approximation of the T operator in the soft dielectric and soft 
diamagnetic Umit as 



(11.4) 



r!=0 



This approximation is more valid for weaker V, so here is 
where the soft dielectric and soft diamagnetic limit is applied. 

Now we will study the lowest expansion order. In the lowest 
expansion order (p — 1 and T — V) we get the asymptotic 
approximation of the Casimir energy between two bodies as 

he f°° 

E = -—j dfcTr(1/iGoi2V^2Go2i). (II.5) 

We can separate the magnetic and electric part of each body 
potential (which is a diagonal operator in positions space, 
because it is local) and each field contribution in the free 
dyadic Green function. For isotropic dielectrics with cons- 
tant electric and magnetic permeabilities, it is easy to find 
that V^{r) = (e„ - eo)Xa(r) = ?aXa(r) and {r) = 
{fJ-a - t^o)Xa{r) = f^aXa{r) in the vacuum, where Xa{r) is 
the characteristic function of the a - body volume (1 inside 
the body and in the rest of the space). Note that we repre- 
sent the electromagnetic properties of the objects by their po- 
tential energy instead the complement representation by their 
boundary conditions ||T^ 1. It will let us to obtain the PSA as 
an integral over the volume of the bodies as required by the 
PSA. The contributions of the free dyadic Green function are 
obtained in the Appendix in Eq. ( IA.19I I - ( |A.22| l. 

As seen in Eq. ( |A.19t - ( |A.22t of the Appendix, operators 
are defined over three different linear spaces: 1) An EH - spa- 
ce, whose components are the electric and the magnetic field; 
2) over the space coordinates, because we are working with 
a vector and a pseudovector; and 3) over positions. We must 
solve the trace over these three spaces: EH - space, vector 
coordinate space and position space. First we solve the trace 
in the EH - space and we obtain: 

E ^ r dkTriVfG^ivfG^i 

I ^/E/-^EH^rH r<HE , T T H r^H E T T E r^E H 

+Vl'G^,WG^,n (11.6) 
We identify each term with this obvious notation: 

E = Eee + Eeh + Ehe + Ehh ■ (11.7) 



A. Purely electric and purely magnetic Casimir energy 

The purely electric case has been solved by K. A. Milton 
et. al. in ifTsll . where they give the following result: 



Eee = — 23eie2 
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(11.8) 



We obtain the same result because the dyadic Green fun- 
ction used in flS*] is the pure electric part of the Green fun- 
ction matrix we use here. In fact, the pure magnetic part of 
the Green function matrix is equal to the pure electric part 
G^j^ {R, k) = Gf;^ {R, k), so the calculus of the purely mag- 
netic part of the Casimir energy is similar to the electric one. 
Consequently we obtain the following result for the purely 
magnetic part of the Casimir energy: 



Ehh = -23/ii/t2 
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To our knowledge, this is the first place where this result is 
obtained. 



B. Coupled electromagnetic Casimir energy 

Here we are going to calculate the contribution of the 
electromagnetic coupling part of the Casimir energy. Using 
Eq. ( IA.2OI 1 and Eq. ( |II.7| i, we have to solve: 

Eeh = r dkTr {V.^G^.WG^i) . (11.10) 
271- 7o 

Replacing each potential by its value and assuming isotropy, 
it is followed that they are proportional to the identity matrix. 
Then we can drop them for the coordinates trace, but not for 
the spatial positions trace: 
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Using Eq. ( IA.20I) of the Appendix and R — |r — r'| we have, 
in matrix form: 

Gi^{R,k) = -k^e^.kdkR 

= Go(i?,fc)A (^k + ^ 

i?3 -R2 
X ( -i?3 Ri \ , (11.12) 

i?2 -Ri 

so the trace in spatial coordinates is easily solved: 

Tr(G^ifG?2f) =-2fc'G2(i?,fc) . (1113) 
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We can simplify our calculus using the dimensionless variable 
u — kR to factorize spatial and frequency contributions of the 
formula 

Eeh = eill2-^a J dri J ^2^, (11.14) 



|ri — 1*21 = R, where R is assume to be a constant. Then the 
integral can be easily solved in this limit to: 
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where 



(11.15) 



Finally, using Jp°° duu^e = , the final result is obtai- 
ned as 
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EH 
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There is an antisymmetry between G^^ and G^^ shown 
in Eq. ( IA.20I ) and Eq. ( IA.21I ) of the Appendix, that is 
G§"{R,k) = ~G"f{R,k). Therefore, we can obtain the 
coupling between the magnetic part of the first object and the 
electric part of the second one in a similar way. Then, we can 
obtain the coupling between the magnetic part of the first ob- 
ject and the electric part of the second one: 



E 



HE 



he 



dridr2 



(4^)3 J J, In-r^r- 



(11.17) 



These two cross contributions to the Casimir energy in the 
soft dielectric and diamagnetic limits are the main result of 
this article. It is interesting to note the sign change in this part 
of the Casimir energy with respect to the other contributions. 
Then, these terms can invert the typical attractive nature of 
the Casimir energy to repulsive in this limit for objects with 
very different electromagnetic nature. Finally, thanks to these 
new two terms, we obtain the complete electromagnetic Casi- 
mir energy between two objects in the soft dielectric and dia- 
magnetic limit for the first time. Then the global asymptotic 
electromagnetic Casimir energy is finally obtained: 
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dridr2 



(11.18) 



where 7 = 23eie2 — 7eifi2 — 7^21^1 + 23/ii/i2 is the comple- 
te multiplicative constant, taking into account all the electro- 
magnetic effects in this limit, nor just electric effects, where 7 
would be 23eie2. 



where Va is the volume of the a - object. In the soft limit order 
we can approximate the electric and magnetic polarizabilities 
as = and = Ut-- In 1 14] it was derived for 
an sphere. By using the same method with the approximation 
of that the effective field over the dielectric is equal to the 
induced field in the soft limit, we arrive at — for any 
arbitrary shaped object. Equation dll.lSb simphfies into 
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rj H E 

"oil 0.2 



(11.20) 



This is the soft response limit of the Feinberg and Sucher po- 
tential 1 10]. It coincide as well with the soft response limit of 
the asymptotic Casimir energy obtained in ]7] for spheres. The 
presented scheme looks as we would use a local two points po- 
tential (where we substitute the polarizabilities of the bodies 
by their local susceptibilities) as in the present case with the 
potential given in Eq. dll.lSI l, and integrate over the volume 
of each body to obtain the PSA of the Casimir energy. If we 
also study the asymptotic distance limit of the PSA, then we 
reobtain the Feinberg and Sucher potential in the soft respon- 
se limit, now proportional to the polarizabilities of the objects 
instead the susceptibilities. 



III. DILUTED LIMIT AT ANY TEMPERATURE T 

In this section we will calculate the Casimir energy in the 
diluted limit at any finite temperature. Then we will focus on 
different approximations to low an high temperature limits, re- 
covering the zero temperature case showed before and giving 
new formulas of the PSA for any temperature. We begin this 
study with the Casimir energy formula jlLlI) for any tempera- 
ture: 



^ = ^5!^^ log(|I-N(fc„)l), (III-l) 

n=0 



C. Asymptotic Casimir energy 



We can also obtain from Eq. ( III.I8I 1 (always in the soft elec- 
tric and magnetic limit) the first asymptotic energy order for 
two objects. These formulas are the first order of a multipolar 
expansion of the integrand in Eq. ( III. 181 ) in the coordinate sys- 
tem of each object. That means that we have to assume that 
the distance between the objects is much greater than their 
characteristic lengths R^, that is R^ <C R- In this limit we 
can separate the problem into two scales and we can replace 



with Matsubara frequencies fc„ = 2n!^^n = An. The til- 
de means that the n — case is weighted by a 1/2 factor 
As usual, we apply the Born approximation to the T ma- 
trix scattering obtaining at first order Tq, — Va (where a la- 
bels the body in interaction). We apply again that log(|yl|) = 



Tr{\og{A)) andlog(l^ 
Eq. ( HITT] ) into: 



-x) 



El — , then we transform 
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E = -kBTj2-J2'^^ (Nf (An)) 



(111.2) 
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At first order in p we obtain: 

oo ' 

E = -kgTY^ Tr(N(An)) 

oo ' 

= -fcBr^Tr(TiUi2T2U2i). 



(III.3) 



As we are working in the positions representation space ins- 
tead in the multipolar representation space, we have to repre- 
sents the operators Tq and Vap in the positions space. For that 
issue, we have into account that Vap matrices represent the 
field propagation between two points. Then they are represen- 
ted by the free vacuum Green function of the interaction field. 
Using the Born approximation, we obtain the next formula for 
the diluted approximation of the Casimir energy: 



-kBTyTriViGoi2V2Go 



(111.4) 



Where R = — r^. The trace in spatial coordinates is easily 
solved obtaining: 



E 



EE 



= -ksT^ eih I dri / dra 



n=0 



,-2kR 



n J2 ^AnYR^ 

X [6 + l2kR + IQk^R^ + 4k^R^ + 2k^R'^] . 

(III.8) 

Replacing k by An we obtain: 



Eee ~ — 



keT , , 

(47r)2'i'2 A,/, Re 



dridr2 ^ ^-2ARn 



n=0 



x[6 + 12Ai?n4- lOA^i^^n^ 
+'iA^R^n^ + 2A^R^n^]. 



Having into account that X]^o ^ 
deduced: 



(III.9) 

it is easily 



di J2 = = ^^(1 " e-^yK (IILIO) 



As we did with ( 111.5b . we can separate the magnetic and elec- 
tric part of each body potential and each field contribution in 
the free dyadic Green function. We trace over the EH space 
obtaining: 

oo ' 

I TrE^EH^rH r<HE I T,rH r<HET,rEr<EH 
+ ^1 ^012 ^2 ^Q21 + ^1 ^012 ^2 '^021 

+V,"Gg^V,"G^,?). (III.5) 
As before, we identify each term with this obvious notation: 
Et = Eee + Eeh + Ehe + Ehh ■ (III.6) 

A. Purely electric and purely magnetic energy 



Using Eq. ( |A.19t and Eq. ( |A.22| i of the Appendix, we can 
solve the trace over the coordinates of Eee and, similarly, 
of Ehh- The matricial form of the purely electric part of the 
dyadic Green function and of the purely magnetic part of the 
dyadic Green function are: 




Rx Ry Rx Rz 

Rl 



RyRz 



R, 



X (i + ZkR+k^R^ 




2u2\ 



X (1 + kR+k^R 
GgfiR,k). 



47ri?5 
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(III.7) 



Denoting A = i?A, we can carry out the sum obtaining: 



E 



EE 
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(47r)' 



reie2 



dridr2 
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i?6 (e2> - 1) 



X [e^^ (-9 - 12A + lOA^ - AX^ + 2X'^) 

+e'^^ (6 + 36A - lOA^ - 12A'^ + 22A4) 

+e^^ (6 - 36A - lOA^ + 12A^ + 22A'^) 

+e*^ (-9 + 12A + lOA^ + AX^ + 2A'*) 

+3ei°^ + 3]. (III. 11) 

We obtain a similar result for the purely magnetic Casimir 
energy replacing 6162 by /ii//2. 

B. Coupled magnetic - electric Casimir energy terms 



We perform the same calculations using (IA.20b and ( IA.21b 
of the Appendix as for the purely electric case. The trace over 
spatial coordinates is already done in (III. 13b , so we obtain the 
formula for the coupling terms of the Casimir energy just re- 
placing the integral in (III. 14b by a sum: 



P 2kBT^ . 



dridr2 



1J2 

2R^k^ + R^k^ 



00 ' 

n=0 



X {R'^k^ + 2R'k-' + R^k'^) . (III.12) 

Here fc„ — 27r-^|^n = An are the Matsubara frequencies. 
After solving this sum and denoting A = RA, we obtain: 

^ ^ dridra A^ 



E 



EH 



2kBT^ . 

T2^lA^2 



(4^)^^^"^7iy2 i?« (e 
x[e2^(l-2A + A2) 

+e-*^ (-1 -6A+ llA^) 

+e^^ (-1 + 6A+11A2) 

+e«^(l + 2A + A2)]. 
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(III.13) 
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Replacing ei/i2 by jlie2, we obtain the formula for Ehe- The- 
se formulas for the Casimir energy in the diluted limit are va- 
lid for any temperature, but they are too much complicate for 
analytical analysis at any temperature. Therefore, we study the 
limits at high and low temperatures. 



C. Low and zero temperature limit 



Where AiE''^' is the first correction to the classic limit of the 
Casimir energy E'^K With these results we see that we loose 
the coupling between electric and magnetic energy terms in 
the classical limit. 

Finally we must remark that, if we take the asymptotic dis- 
tance approximation as in Sect. Ill CI we reobtain the results 
given in [Ts^ in the soft response limit. 



Here we will recover the zero temperature limit of the di- 
luted limit. Then we will make a perturbative analysis valid 
for low temperatures. We can make it easily by expanding the 
Taylor series in A = RA of the integrand and studying just 
the first orders deleting the rest ones. By using dlll.l 11 ) and 
( IIII.13I ). we reobtain ( III.I8I 1 for the zero temperature case. We 
need to take the fifth order series term of the Taylor expan- 
sion of Eq. dlll.l 11 1 and Eq. ( IIII.13I ) to get the next non zero 
perturbation term of the Casimir energy as 



£162 



227r3 fksT 
-kbT 



945 



he 



dridr2 

1J2 R ' 
(III. 14) 



dYldY2 



A^Eeh _ ^5Ehe _ 27r3^^^^ /fc^r 

eiA2 h^2 189 \ fic J J1J2 R 

(III. 15) 

because AiE = A2E = A3E = A4E = 0, where AnE is 
the n-th order correction term to the low temperature expan- 
sion. 



IV. PSA FOR THREE BODIES SYSTEM 

In this section we are going to calculate the Casimir energy 
between three bodies in the Pairwise Summation approxima- 
tion. In this asymptotic limit we will obtain that the energy 
of the system will be the addition of the PSA energy of each 
pair of objects. This linear behavior of the Casimir energy was 
expected in that approximation although we know that it is in 
general false. We begin this study from the Casimir Energy 
formula for three objects given in [16.1 : 



E=^ I dfclop 



(IV. 1) 



The AI matrix (whose coefficients are non commutative ma- 
trices) is: 



^-1 


-U12 


-Ui3 


-U21 




-U23 


-U31 


-U32 





D. Higli temperature and classical limit 

In this section we will obtain the high temperature limit of 
the Casimir energy in the diluted limit, whose first term will be 
the classical limit of the Casimir energy. To obtain this limit, 
instead of solving the sum in Eq. (Illl.St and Eq. ( |III.12| i, we 
will just keep the first term of the sum. Then the classical limit 
of the Casimir energy in the diluted limit reads the first sum 
term: 



Eee 
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fllfl2 



' (47r)2 



dridr2 
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(III. 16) 



(III. 17) 



ei/i2 Mie2 

And the first perturbation to that limit is the next sum term: 

^lE^E ^ _ML f f dridr2 ^_2KR 



WJJ2 



(6 + 12Ai? + lOA^i?^ + 4A^i?^ + 2A^i?^) 

•cl 

(III. 18) 
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fiie2 (47r)2 
X {R'^A'^ + 2R^A^ 



dridr 



■1J2 R' 
R^A^) . 
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(III. 19) 



So, using the logarithm product rule, the Casimir energy bet- 
ween three bodies is 

he r°° 

E3 ^ — dfcl0g(|I-Ni2|) 

27r Jo 

) 

&l0g(|I-Ni3|) 



he 
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+ ^^°°dHog(|I-R|). 
Here the A'' matrix is Nap = 'T^aVap'^iiVpa, and: 



(IV.2) 



M = (I-Ni3)"'(T3U32+N3i) 

X (I-Ni2)"'(T2U23+N2l). (IV.3) 

Expanding the inverses we obtain the double series 

00 

M = ^ (Ni3)"^ (T3U32 + N31) 

ni=0 
00 

X (^12)"^ (T2U23 + N21) . (IV.4) 

"2=0 

The lowest order of the Born series of the M matrix will come 
from the first order expansion series making ni — n2 = 0, 
that is 

M~ (T3U32+N3l)(T2U23 + N2l). (IV.5) 
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In addition to that, we just consider the sum term with the 
minimum number of T matrices products, because that will be 
the lowest order expansion in susceptibilities of the M matrix. 
That means that we reduce the highly non linear M matrix to 

K~(T3U32T2U23) =N23. (IV.6) 

Replacing Eq. ( IIV.6I 1 in Eq. ( IIV.2l i. we obtain the next PSA of 
the Casimir energy between three objects: 



he r°° 
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dfcl0g(|I-Ni3|) 



he r°° 

+ — *l0g(|I-N23|). (IV.7) 

Jo 

Therefore we obtain that the PSA approximation of the Casi- 
mir energy between three objects is the sum of the PSA energy 
of each pair of objects. In other words, we rederive a kind of 
superposition law of energies in the diluted PSA limit as ex- 
pected, because the usual presentation of the PSA approxima- 
tion is the assumption that we have a superposition behavior 
in that asymptotic limit. Here we have justified the validity of 
that approximation. 

If we take the asymptotic distance limit to Eq. dlV.TI l, we 
will obtain a superposition of two bodies PSA energies in the 
asymptotic limit. The nonlinearity of the Casimir energy must 
be given by higher orders expansion terms even for systems 
with three bodies. The same superposition behavior is expec- 
ted for the general N body case, and it will be proven in the 
next section. 



A. PSA for general N bodies system 

In this section we are going to generalize the PSA energy 
of three bodies given by Eq. (IIV.7l i. For this purpose we are 
going to use an iterative procedure which will give us the new 
terms to include to the PSA of the Casimir energy of n — 1 
bodies when we include another new object to our system. 
Let us represent the M matrix of Eq. dlV. ll) as the sum of its 
diagonal and its non-diagonal parts as iflql : 



or symbolically as 



■U, 



(IV.8) 



(IV.9) 



then it is easy to find that the inverse of M is the next pertur- 
bative series 



-1 = (-UT)' 



(IV. 10) 



)i=0 



On the other hand, the M matrix of the iV objects system is 
related with the M matrix of the — 1 objects system by block 
matrices in the next way 



-U^,Ar-l TT^ 



Where 7 index goes from 1 to iV — 1. With this result we will 
calculate the determinant of M^v obtaining 

IMatI = |MAr_i||TAr|-l|I-TjvU^,Ar_iM^l_^UAr_i,^|. 

(IV. 11) 

In the PSA approximation is valid the substitution M^^ ~ T, 
where T is a diagonal matrix whose — 1 diagonal elements 
are T-^ with 7 index defined as before. Then, using this appro- 
ximation we can approximate the Ma? determinant as: 

IMatI - |MAr_i||TArrl|I-T^U^,Ar_iTUAr_i.^|. (IV12) 

Where we multiply by blocks the matrix of the last determi- 
nant obtaining 



IT 



/NiTv 

N2Ar 

V 



Finally, we obtain the desired result: 



a AT 



N\ 



7=1 



The other needed matrix is 



IT 



wl 



And with the initial conditions: 

|Mi| = |Moo,i| = ITTil 
we obtain these determinants in a closed form as 



1-1 



(IV. 13) 



(IV. 14) 



(IV 15) 



(IV 16) 



(IV. 17) 



fe=i 



and 



i=2 m=l 



AT 



(IV. 18) 



Taking the logarithms of Eq. ( IIV.II ). the Casimir energy for the 
N bodies system is approximated as 



JV (-1 



he f°° 

En^— dfc^ ^ log(|I-N,,„|) (IV19) 

in the PSA limit. Eq. ( |IV.19l l is the generalization of Eq. ( IIV.7I ) 
for the N body case, and it show us that in the PSA limit, we 
always obtain a superposition principle of the two body PSA 
Casimir energy despite the fact that Casimir energy is not a 
nonadittive interaction. We will obtain the non-linear effects 
in the next orders of the expansion. 

This result is qualitatively different of the usual PSA proce- 
dure of integrating each point of a point potential over the 
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volume of each object. The reason is simple, a N point poten- 
tial is proportional to N polarizabilities lUTIl . but in the soft 
limit approximation, the lower allowed term is proportional to 
two polarizabilities. Therefore we should study the {N — 1) 
expansion term to obtain the first one proportional to N pola- 
rizabilities. So in a consistent calculation of PSA energies for 
three or more objects, the asymptotic approximation of the 
PSA energy in the soft material approximation is different of 
the N point potential function. In the next section we will find 
where the contribution of the three point potential is relevant 
in the diluted limit. 



V. SECOND ORDER EXPANSION OF DILUTED LIMIT 

In this section we are going to study the second order ex- 
pansion of the PSA in the diluted limit. This is a complica- 
ted long calculus so, instead the complete electromagnetic ca- 
se, we will restrict ourselves to the purely electric case. That 
means that we will ignore the magnetic properties of the bo- 
dies. It could be interesting to show the complete study of that 
series term, but it is a long calculation to show here. However, 
it possesses theoretical utility, as it shows the lack of the su- 
perposition principle in the Casimir energy calculations. The 
next results have been done with the help of Mathematica 
CI. 

Again we start from the Casimir energy between two com- 
pact objects as given by Eq. ( III.2I I. We studied the first order 
expansion of the energy in Sect. |ll] Now we are interested in 
the next order expansion term in polarizabilities, so we main- 
tain the study of the fist term of Eq. ( III. 2b taking just the case 
p — 1. In addition to that, we take the second order approxi- 
mation of the Born series of the T matrix: 



Where k is the frequency, a and /3 labels the bodies coordi- 
nates and i shows the vector component of R^^ used. So we 
can solve the trace over the space coordinates obtaining the 
following result for the first term of ( IV.2I) : 



dri / dv: 



he 

„-k(Ri2+R2l'+Rl'l 



2 / dry 

2 Jl' 



-€16261/ 



dk 



x[7^[l,2],(fc)7^[2,l'],(fc)7^[l^l],(fc) 
+7^[l, 2]y(fc)7^[2, l%(fc)7^[l^ i]y(fc) 
+7^[l, 2],(fc)7^[2, i'Uk)n[i\ lUk)]. 



(V.4) 



Here the space integrations are performed just over the body 
volumes because it is there where dielectric and diamagnetic 
potentials are defined. It is possible to perform the k integral, 
because we have an exponential multiplied by a polynomial 
integrated between zero and infinity. It is not shown here be- 
cause it contains 216 sum terms. 

For any finite temperature we also obtain a long result, but 
in the classical limit we can at least obtain a tractable formula. 
This limit consists of keeping only the fc = mode. In that 
case we obtain the result 



he 



dri / dr2 / dvi 



616261' 



256^4 --j^, Rl^R^^M^,, 

x[7^[l,2],(o)7^[2,l'],(o)7^[l^l],(o) 
+7^[l, 2]y{o)n[2, i']y{o)n[i', i]yio) 

+7^[l, 2],(0)7^[2, l%(0)7^[l^ 1],(0)]. (V.5) 



Where 7^[a,/3],(0) - ( i?^ ^ 



SO we can write: 



n=0 



y ~ F + VGnV. 



(V.l) 



Applying the linearity of the trace and taking our attention just 
to the second order expansion in polarizabilities (that is, to the 
third order term in polarizabilities), we obtain the following 
result for the second order of the diluted limit result: 



he 

he 
'2^ 



— / dkTTiVyGoi'lViU,2V2U21') 



dkTl{ViUi2'V2'G„2'2V2U2l). (V.2) 



We will center our study just on the first integral of ( |V.2| l. be- 
cause the analysis of these two integrals is the same. First we 
make the trace over the EH space, which is automatic here 
because we have cancelled the magnetic properties of the bo- 
dies. In other words, we are not studying the coupling between 
electric and magnetic induced dipoles. After that we trace over 
the space coordinates, so we need the matricial form of the 
electric part of the dyadic Green function given in Eq. ( |A.19l l 
of the Appendix, where R = r — r'. We will need to define the 
function: 



7^[a,/3],(fc) = Rlp{l 



+ kRal3 

3kRaj3 4 



k^R 



r2 



(V.3) 



Et. = - 



he 



2567r4 
x[(i??2- 



dvi / dY2 I dvi 



eie2ei' 



RO PO P5 
-n-12-'T-21'-'^l'l 



3^122:, 



1-12 



3Ri 
3i?i2 



2y, 



{R21 
(Rl 



21' 
2 



3R: 



3R 



21'x 
2 N 



I (i?^i 



- ^Ryix) 
^Rviy) 

^R\'Iz)]- 

(V.6) 



When performing this integral, a term proportional to Ry\ 
and another one to R^i^ appear These terms can be proble- 
matic because i?i'i = belongs to the integration interval 
for all the points of the 1 body. But these singularities can be 
removed with the appropriate regularization procedure. No- 
te that if we make 1' = 3 in Eq. ( |V.2| i. this integral would 
also be one of the contributions of the second term expan- 
sion of a three objects system. In that case we have not got 
any problem in the volumes integration for this term and this 
expansion term looks quite similar to the local analog of the 
three points potential given in |17]. In fact, in the three (or 
more) bodies system that term is the first term that breaks the 
superposition behavior founded in|lVl Replacing the fc^ term 
by the operator - A in Eqs. ( IA.19I ) and ( IA.22I ) of the Appen- 
dix, and replacing the k term by the operator V^A in Eqs. 
(IA.2II 1 and ( IA.20b of the Appendix, which is a change valid 
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on shell because Gq{R, k) = , we recover the forma- 
lism used in Iil7i1 in the soft dielectric and soft diamagnetic 
limit. This scheme let us take into account not only electric 
phenomena as in [ 17], but also magnetic and electro-magnetic 
coupling effects. Taking this argument into account, we can 
understand the structure of the perturbation terms of the PSA 
limit of the electromagnetic Casimir energy between N bo- 
dies. The Casimir energy in the PSA limit has the structure 
of a series of infinity terms whose n-order term is the sum of 
integrals over rt + 2 bodies (which can be repeated or not) of 
the [n + 2) points local em potential. Each series term come 
from the expansions made in Eq. ( III.2I ) and in Eq. ( III.4I ) and 
is proportional to the product of n + 2 permittivities. If we 
take the asymptotic distance limit of these integrals, we will 
obtain a series of n points em potential now with polarizabili- 
ties instead susceptibilities as in [17], but in the soft response 
limit. When we use perturbations of the tree term of the Born 
expansion, it appears divergent terms in the integrand as seen 
in Eq. (IV.6b . These singularities can be removed with the ap- 
propriate regularization procedure. 



VI. FINAL REMARKS 

We have calculated the full electromagnetic Casimir energy 
in the diluted limit between two bodies using the formalism 
given in ]7] reobtaining the energy given in ]6] in the pure 
electric case. Although this formalism fails in the perfect me- 
tal case, it is valid for soft dielectric and diamagnetic bodies. 
This formalism can be extended for an arbitrary number of 
bodies. We have shown that the Casimir energy, which is a no- 
nadittive interaction, has a superposition behavior in the first 
expansion of the PSA and how the integration of different N 
point potentials appear in the perturbation series to contribu- 
te in the whole energy. Then we have a systematic procedure 
to calculate Casimir energies valid for soft dielectric and dia- 
magnetic bodies. 
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Appendix: Obtention of the matricial Green function 

Here we have calculated the Casimir energy using a matri- 
cial Green function for the electromagnetic fields instead the 
dyadic Green function used in ]'7']. in this appendix we will see 
that both formalisms are equivalents just using as field sources 
the induced polarization and magnetization vectors instead the 
induced currents into the bodies. The induced currents into the 
bodies ai^e 1,19,1 : 



We introduce this linear change of variable in the partition 
function, so there is not any new relevant term in the action of 
the problem. This action will be transformed in the next way: 

S = j rixj^oj, (A.2) 

dxfc^P^oP + j dxV X M^oV x M 
dxikPGoS/ X M - y" dxV X MGoikPiA.S) 



S 



Where 



Go 



1 ^ ^ 



Go, 



(A.4) 



is the Green dyadic function. We integrate by parts each action 
term obtaining: 



See 



dxk^FGoF, 



S 



EH 



dxPiCfj^Pj, 



dxikPGo^ X M, 



(A.5) 



/ 



= / dxikPi 



1 



(A.6) 



Seh ^ J dxP^[ike^c.pV^]GoW 

+ J dxP, ifc^VjVjejc^/jV 



GoM".(A.7) 



The second term is zero: 

eja/jVjM^V^Go = M • V X VGo = M • = 0. (A.8) 
So we finally get: 

Seh = J dxP, [-ike^jpV^] GoM^ 

dxPiGf^M'. (A.9) 



She is similar to Seh, but with a changed sign because the 
complex conjugation of the induced current into the action: 



She = - / dxV X MGoikP, 



(A. 10) 



-ikP + W X M. 



(A.l) 



She = j dxM.Gg^P^ = J dxM, [-Gf^"] P^ (A. 11) 
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And the last term: 

Shh = / dxV X M^^oV X M, (A. 12) 



T^e^"''e,a6V«V,V'^ V"'Go - 0, (A. 17) 



Shh = J dxe'"^V„M^ 



Shh = I dxe'^'^M^V 



(A.13) 



(A. 14) 



because this differential operator is zero. So Shh is 



Shh = ydxM4e'"^V^%Goe,abV'^)M'^ 

-J dxM« (^e-'3^^^V,V'^Goe,ahV'^^ M^ 

(A. 15) 

Using e'^^Sfejab = <5?(5f - S^S^ and VhV'^Go = PGq, we 
get the next result; 



= [j^VfcV^- VaV'"] Go, 

= [<5^A;2 - VaV'"] Go. (A. 16) 

The other term requires even a more tedious work, but it is 
easy to obtain that 



5* 



HH — 



y"dxM„ [S^^k^-VaV"] GoM°, 



(A. 18) 



After a Wick rotation, we obtain the used form of the matri- 
cial Green function, which components are, using Go {R,k) = 



AttR 



andi?= |r-r'|: 



Ggf{R,k) 
Ggf{R,k) 



Which are the results used in Eq. (|II.5t . It is also possible to 
obtain the same result from fluctuation - dissipation theorem, 
as made in ll20ll. 



[fc2%+V,V;] Go(i?,fc), (A.19) 

-fce,,fcVfeGo(i?,fc), (A.20) 

fce»jfcVfcGo(i?,fc), (A.21) 

[fc2j,^-+V,V;-] Go(i?,fc). (A.22) 
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